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Abstract 

We present a predictive scheme for fermion masses and mixings in- 
spired by supersymmetric SO(10) in which the gauge hierarchy prob- 
lem is resolved without fine tuning the parameters. There are six 
predictions in the flavor sector, all consistent with the present data. 
The scheme reproduces the familiar asymptotic relations m& = m T 
and mdm s mb = m e m^m T - A new expression for V c b is obtained in 
terms of the quark masses. The remaining predictions involve the 
quark mixing angles V us and V u b, as well as the parameter tan/3 which 
turns out to be close to mt/mb. 
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1. Introduction 



In this paper we present a scheme for fermion masses and mixings motivated 
by supersymmetric SO (10) grand unification (GUT) in which the gauge hier- 
archy problem is overcome without fine tuning. There are six predictions in 
the flavor sector, all of which are phenomenologically consistent and compat- 
ible with the existence of a heavy top quark. A simple derivation in SO (10) 
of the quark and lepton mass matrices that lead to these predictions will be 
presented. 

The familiar asymptotic relations = m° T and m°m°m° = m° e m ^m° T 
(see [l])g are reproduced in our scheme, while a new sum rule for m° in 
terms of the quark mixing angle V c b is obtained. Two of the remaining three 
predictions are for V us and V u b, resembling those from the Fritzsch scheme [2], 
while the third prediction is for the ratio of the two Higgs vacuum expectation 
values, tan/3 ~ mt/rrib [3]. 

The model presented here provides a unified framework for some of the 
well-known and successful relations for the quark masses and mixing an- 
gles which are compatible with a heavy top quark. Since the up-quark, 
down-quark and the charged-lepton mass matrices are assumed to have sim- 
ilar forms, their derivation from an underlying SO (10) GUT is relatively 
straightforward. 

The SO (10) model is the simplest grand unified example that assembles 
all fermions of a given family into a single irreducible representation, the 
spinorial 16. As such, it has proven useful to work with SO (10) in attempting 
to model fermion masses and mixings. It also has the potential to relate the 
neutrino masses and mixings to those of the quarks and charged leptons 
[4]. The doublet-triplet mass-splitting problem that all GUT models must 
address has an elegant resolution in SUSY SO (10) without any fine tuning 
of the parameters [5,6,7]. We are also encouraged by recent attempts to 
derive realistic supersymmetric SO (10) GUTs from superstrings in the free 

1 We use the superscript to denote GUT scale quantities. 
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fermionic formulation [8,9]. 

In the 5*0(10) model of fermion masses that we consider, we shall re- 
quire that the doublet-triplet mass splitting problem be resolved without 
fine tuning. Furthermore, motivated by our desire to preserve the successful 
prediction of sin 2 ^, we assume that SO (10) breaks directly to the minimal 
super symmetric standard model (MSSM). These two requirements typically 
call for a non-minimal Higgs system at the GUT scale which we exploit to 
arrive at a predictive scheme for quark and lepton masses and mixings. 

2. Model 

We shall adopt the mechanism developed in Ref. [6,7] to solve the doublet- 
triplet splitting problem naturally in SUSY 5*0(10), making use of an old 
suggestion of Dimopoulos and Wilczek [5] . This involves the coupling of Higgs 
10's with an adjoint 45: 10i. 45.102- If the 45 acquires a vacuum expectation 
value (VEV) along the (B — L) direction, (45) = diag.(a, a, a, 0, 0) x ir2, 
this coupling gives GUT scale masses to the color-triplets in the 10's, while 
leaving massless the SU(2) L doublets to be identified with H u and H d of 
MSSM. The superpotential term (MgutIO 2 .) will make the doublets in 10 2 
superheavy, leaving MSSM as the effective low energy theory. 

The minimal Higgs sector that can break the 5O(10) gauge symmetry 
down to MSSM involves one 45 and a spinorial 16 + 16 [7]. However, the 
low energy theory below the GUT scale in this case is not quite MSSM. In 
particular, some of the Higgs(ino) superfields turn out to have intermediate 
scale masses. Here, since we require that the theory below the GUT scale be 
MSSM, it will be necessary to have a slightly more extended Higgs system 
at the GUT scale. 

The Higgs system that we employ for the symmetry breaking involves two 
sectors. One sector breaks 5O(10) down to SU(5). This can be achieved for 
example, via the following superpotential involving a 45 (A) and a 16 + 16 
(C + C): 

W 1 = m 1 CC + m 2 A 2 + X 1 CCA . (1) 
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This superpotential induces a VEV for C and C along the SU(5) singlet 
direction, and for A along SU(5) x U(l)x' (A) = diag.(a,a,a,a,a,) x ir 2 . 

The second sector breaks 50(10) down to SU(3) C x SU(2) L x SU(2) R x 
U(1)b-l- The superpotential involves a 54 (S) and a second adjoint A' [10]: 

= m 3> S 2 + m 4 A' 2 + A 2 S 3 + X 3 A' 2 S . (2) 

The resulting VEVs of S and A' are (S) = diag.(s, s, s, — §s, — |s) x 1, (A') = 
diag.(a', a', a', 0, 0) x ir 2 . The A' field in this sector is also responsible for 
the doublet-triplet splitting via the superpotential terms involving Higgs 
10-plets (denoted by H and H') 

W 3 = HA'H' + m 5 H' 2 . (3) 

In order to keep the masses of the electroweak Higgs doublets at the 
correct scale, the VEV of A' should be along (B — L) to a high degree 
of accuracy. This requires that the A' coupling to the (A, C) sector be very 
weak, since any such coupling will induce an A' VEV along the SU (5)xU(l) x 
singlet direction as well. Now the (S, A') sector should be linked with the 
(A, C) sector, otherwise pseudo-Goldstone bosons will result. The simplest 
way to link the two sectors without upsetting the VEV pattern of A' is [6] 
by a coupling involving another adjoint: Tr(AA'A"). This term, due to its 
complete antisymmetry, vanishes at the minimum, and thus does not affect 
the VEV of A'. Yet, it gives GUT scale masses to all the would-be pseudo- 
Goldstone bosons. The A" superfield can have its own superpotential and 
its VEV can in general be written as (A") = diag.a"(l + z, 1 + z, 1 + z, 1 — 
|z, 1 — |z) x ir 2 . Here z = would correspond to the VEV being along 
U(l)x, while z — — 1 (2/3) would preserve I 3 r (B — L). 

Since our focus here is the fermion Yukawa sector, we shall not get into 
the details of symmetry breaking, except to note that one can find discrete 
symmetries that would lead to the superpotential given in Eqs. (l)-(3) along 
with the Ty(AA'A") term, while preventing other terms which potentially 
can upset the gauge hierarchy. Note that one has the option of using gauge 
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singlet superfields with GUT scale VEVs to induce some of the mass terms 
in Eqs. (l)-(3). 

We now proceed to show that the above Higgs structure can lead to a 
predictive set of mass matrices for the quarks and charged leptons, provided 
one introduces a flavor-dependent discrete symmetry. 

3. Fermion Mass Matrices 

In our scheme, only the third family receives mass from a renormalizable op- 
erator in the superpotential. The second family masses as well as the mixing 
angle arise from dimension 4 operators in the superpotential, suppressed by 
one inverse power of M, where M is a scale much larger than the GUT scale. 
The first family masses and mixings will come from dimension 5 operators. 
It is entirely conceivable that these non-renormalizable operators arise from 
integrating out some vector-like families with masses of order M. The inter- 
generational mass hierarchy is thus related to the small ratio (Mqvt/M), an 
idea discussed by numerous authors in the past [11]. 

In our discussions we shall not assume any special structure for these 
'effective' operators. They will only be constrained by the SO (10) gauge 
symmetry and some discrete flavor symmetry. In particular, if one contrac- 
tion of the group indices is allowed in a given nonrenormalizable operator, we 
shall allow for all possible group contractions. This should make the deriva- 
tion of such 'effective' operators from an underlying theory somewhat easier, 
since no special care is needed in the way the heavy fields are integrated out. 

As indicated above, the third family acquires its mass from a direct cou- 
pling to the Higgs 10-plet H . All the other masses will involve the adjoint 
fields A and A" with VEVs specified earlier. In order to generate the first 
family masses and mixings, a singlet field X with a GUT scale VEV is also 
employed. The Yukawa superpotential is given by 

WVuk = h^C la ^H a + h 2Z i) T 2 C 1{ahc] i) Z A" ah E c + h' 2 ^lC la ip z A^ b H b 
+ h 12 ^C 1{abc} ^ 2 A ab H c X + h' 12 ^C^ 2 A ab H b X . (4) 
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Here fy, (i — 1 — 3) stand for the three fermion families belonging to the 
16 of 5*0(10), C is the SO (10) charge conjugation matrix, and 7{ a fe c } is the 
totally antisymmetric combination of the SO(10) gamma matrices: 

l{abc} = lalblc ~ lalclb ~ Iblalc + Iclalb + Iblcla - Iclbla • (5) 

In Eq. (4), appropriate powers of M in the denominators are to be under- 
stood. 

As an example of a discrete symmetry that would lead to the Yukawa 
terms of Eq. (4), consider the transformations of the relevant fields under z 5 : 
(ipi, fps) ~ (u 2 ,u,l), A" ~ uj a ,X ~ lu 2 , where u 5 = 1. This guarantees 
the radiative stability of the Yukawa couplings. Note that in the higher 
dimensional operators in Eq. (4), both the 10 and 120 'effective' operators 
contribute, i.e., all possible contraction of the SO (10) group indices have 
been included. 

The couplings h 33 , h 23 , hi 2 can be made real by field redefinitions, while 
h' 12 , h 23 and z are in general complex parameters. 

To arrive at the charged fermion mass matrices from Eq. (4), it is neces- 
sary to determine the couplings of all the components of 16j with the light 
doublets (H u , Hd) in H. This is carried out by adopting a definite set of 
SO (10) gamma matrices [12]. The fields H u and Hd in the notation of Ref. 
[12] are H u = (H 9 + iH 10 )/V2 and H d = (H g -iH 10 )/y/2. Consequently, the 
/123 term in Eq. (4) can be expanded as (after absorbing a" /M into /123, and 
similarly (AX) /M 2 into h 12 ) 





V2ih 23 H u 


u 2 u c 3 (- 


~Qu + 2Qu c — (1 - 


3 

- -z))+u 3 u c 2 ( 


—Qu c + 2Qu — (1 


- 2 Z)) 


+ 


V2ih 2 zH u 


v 2 vl(- 


-Q u + 2Q u c - (1 - 


3 

- -z)) + u 3 u c 2 (- 


-Q„c + 2Q U - (1 - 


— z)) 
2 " 


+ 


\f2ih 2 2,H d 


d 2 dl(- 


-Qd + ZQdc + (1 - 


lz)) + d 3 d c 2 (- 


-Qv + 2Q d + (1 - 


- -z)) 
2 JJ 


+ 


V2ih 23 H d 


e 2 e c 3 (- 


-Q e + 2Q e c + (1 - 


3 

-z)) + e 3 e c 2 (- 


Q eC + 2Q e + (1 - 


3 \<\ 

2 Z) \ 



(6) 
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Similarly, the h 23 term gives rise to 
3 

V2iti 23 (l--z) [~H u (u 2 u c 3 + u 3 u c 2 + v 2 v c 3 + v 3 v 2 ) + H d (d 2 d 3 + d 3 d 2 + e 2 e c 3 + e 3 e 2 )} . 

(7) 

Here we have defined the "charges" of quarks and leptons to be Q = X + 
6 (j2~), where X is the U(l)x charge normalized such that the 16 of SO (10) 
decomposes into SU(5) x U(l) x as 16 — > 1 5 + 5_ 3 + 10i. These charges are 
given by 

Q u = 1 + z, Q u c — 1 — Az, Q d = 1 + z, Q d c = -3 + 2z, 

Q e = —3 — 3z, Q e c — 1 + 6z, Q u = — 3 — 3z, Q u c = 5 . (8) 

We can now write down the mass matrices for the up-quark, down quark 
and charged leptons. Let us define a 23 = y/2ih 23 , r 23 a 23 = y/2ih' 23 (l — 
a±2 = \/2ihi 2 ,ri 2 = \[2%h! X2 . The mass matrices are then given by 

/ ~a 12 r 12 \ 

M u = v u \ -a 12 r 12 -023(^^ + ^23) 

\ -a 23 (-^z + r 23 ) a 33 J 



M d = v d 



( ° 

012(6 + r 12 ) 

V 



oi 2 (-6 + r 12 ) 


023(6- §2 + r 23 ) 





3, 
2' 

033 



\ 



o 2 3(-6 + \z + r 23 ) 



J 



( o 12 (6 + r 12 ) \ 

M l = v d \a 12 (-6 + r 12 ) a 23 (6 + f z + r 23 ) (9) 

\ a 23 (-6 - y z + r 2s) o 33 / 

A few remarks about (9) are in order. By construction, they have a 
Fritzsch-like texture. However, there is an important difference in that the 
matrices are not symmetric. This is because both the symmetric 10 and the 
antisymmetric 120 'effective' operators contribute to the (12) and (23) entries 
in (9). The net contribution is then neither symmetric nor antisymmetric. 
This implies in particular that the top-quark mass in our present scheme can 
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be larger than the 150 GeV limit [13] set by the Fritzsch ansatz (the limit 
is < 90 GeV without the renormalization group considerations [14]). The 
matrices in (9) admit a top quark as heavy as about 200 GeV. 

Another point to note is that the VEV of A", being proportional to 
the parameter z, breaks the SU(5) symmetry. Thus the bad predictions 
of minimal SU(5), namely m° = mJJ and mj} = m° e will be corrected. 

As noted earlier, the parameters (z, r i2 , r 23 ) are in general complex. In 
the analysis that follows, we shall assume that CP is a good symmetry of the 
Lagrangian. It can be spontaneously broken by the VEV's of A, A" or X. It 
is easy to see that if CP is a good symmetry of the superpotential in (l)-(3), 
it admits a solution where (A) and {A") are real. If this solution is chosen, 
the only possible source of CP violation in the mass matrix is in the phase 
of the singlet field X. If two such fields are present, there is, in general, a 
phase in r 12 which cannot be rotated away. We shall assume in the analysis 
that follows that r 12 is complex, while (z, r 23 ) are real. 

The assumption of spontaneous CP violation is of course motivated by the 
desire to reduce the number of arbitrary parameters. Within the context of 
supersymmetry, there is another reason. It is a well known problem that the 
new phases in the soft SUSY breaking sector of the MSSM, unless somehow 
suppressed, will lead to unacceptably large values of the neutron and electron 
electric dipole moments. One way to suppress these potentially dangerous 
contributions is to assume that CP violation has a spontaneous origin [15]. 

With only r 12 complex, there are 8 parameters (a 33 , a 23 , a i2 , r 23 , r 12 , a, z, v u /vd) 
in Eq. (9) (a is the phase of r 12 ) to fit 14 observables (6 quark masses, 3 
lepton masses, 3 quark mixing angles, one CP phase and tan/5). This results 
in 6 predictions in the flavor sector as advocated. 

Two of the predictions of the model are 

m° b ~ m° T ; tan/5 ~ — (10) 

which can be seen by considering the (33) elements of M u>d j. The resulting 
6-quark mass at low energies is known to be consistent with data [3]. 
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Another prediction of the model is 



(11) 



which follows from the determinants of and M/. This is one of the Georgi- 
Jarlskog relations [1], which also fits the low energy data quite well. 

The remaining predictions have to do with the three quark mixing angles. 
The first one relates V c b to ratios of quark and lepton masses. This relation is 
new, and so we shall explain its derivation in some detail. To arrive at it one 
can safely ignore the masses and mixings involving the first family. Define 



P = Vp 



Q = VQ 



0™0 



mlm 



(12) 



where r] PtQ 
see that 



±1, depending on the signs of the fermion masses. One can 
2(1 - Q) 



r 23 ~~ 



3 + 7Q - 10PQ ' 
225 [(1 + Q- 4PQ) 2 - 4Q] 
(3 + 7Q - 10PQ) 2 

Numerically, the parameters \P\ and \Q\ lie in the range 

\P\ = (— to —); \Q\ = (—to — 
1 1 V5.4 7.0/ 11 V2.9 3.7 



(13) 
(14) 



(15) 



where we have used the relation m° s m® = mW, extrapolated the light quark 
masses from low energies to the top mass scale (~ 160 GeV — 190 GeV) using 
a s (M z ) = 0.12, and imposed the phenomenological constraint m s /m d = (15 
to 25). 

Since r| 3 > 0, 

\1 + Q-4PQ\ >2y h gj. (16) 
Depending on the sign factors r\ P and t]q there are then four possibilities, 
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(i) For r/p — +,T]q — +, Eq. (16) translates into an inequality^ 



> 



41^1 



1 - 



1/4 



2 ' 



(17) 



This leads to a lower limit on m\ hys of about 165 GeV. 
(ii) If r\p = -, T)q = -, then 



\P\ < 



(1-\Q\- 2j\Q\) 



(18) 



which cannot be satisfied since the right hand side is negative, 
(iii) If r]p = +,rj Q = -, 



\1-\Q\+4\PQ\\>2J\Q\ 



(19) 



leading to the constraint 



< 



4 


T° 






fl°d° 


1/4 + 


11° d° 


1/2 - 1) 



(20) 



This results in an upper limit on hys of about 125 GeV which is inconsistent 
with the data. 

(iv) Ifrjp = -,t)q= +, Eq. (16) leads to the inequality 1 + \Q\+A\PQ\ > 
2 y\Q\ , which is automatically satisfied. 

The allowed solutions therefore are r\ P — ±, t]q — +. Using these sign 
factors, the asymptotic expression for V c b can be readily obtained. It is given 
by the relation |V^| = |a 23 /a33||6— 9z-2r 23 |, with \a 2 3/a 33 \ = |c°/t°| 1/2 /|r| 3 - 
(152/2) 2 ! 1 / 2 . Using (13), (14), this yields (x = |e°s%u c/°| 1/2 , y = \c T /t°fjP\) 



1 

V2 



ft 



1/2 



2x — 2r]py ± ((1 + x — Arjpy) 2 — Ax) 



1/2 



(1 + x -2r lP y) 1 / 2 



[21] 



2 For brevity, we denote m® by t and so on. 
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Here we have used the relation \m° s /m Q 



I 1 / 2 which is a consequence 



of the identity |m^m°' 



|. The ± sign in (21) corresponds to choosing 
r 23 to be =)=. Only the positive sign for r 2 3 will lead to an acceptable V^. 

It is instructive to approximate (21) in the limit of an infinite top-quark, 
even though the finite top mass effects will turn out to be significant. In this 
limit, \V® b \ can be written as (for r 2 3 positive) 



\V° h \ ~ — 



1/2 



y/TTx 



(3x 



Vp 



c°r° 



1 + I2x + 3x 2 
1 — x 2 



(22) 



Note that in the strict Georgi-Jarlskog limit, viz., — 3|s°|, = 3|e°|, 
x = 1/3 and the dominant (3a; — 1) term in (22) vanishes. As a consequence, 
the model admits very low values of \V c b\. 

Another useful expression for \V® b \ which resembles the Fritzsch relation 

is 

t° 



[7 


6- 


\z - r 23 


1/2 


w 


6- 


\z + r 23 





f z + r 23 



1 ^ 

4z - r 23 



1/2 



(23) 



It is evident from (23) that for r 23 > 0, the ^Js°/b° term has a suppression 

factor while the ^Jc°/t° term is enhanced (z ~ +1/4), thereby yielding values 
of jV^ | that are significantly smaller than those from the Fritzsch ansatz. 

Relation (21) for \V^\ can be extrapolated to low energies by using the 
renormalization group equations corresponding to tan/? ~ m t jm b . The run- 
ning factors for the relevant quantities to go from the weak to the GUT scale 
can be expressed analytically as (tan/3 ~ m t /m h ) 



Vkm 



Vh/t = 




; Vs/b = Vc/t 



1879/37224 - n2/3 

olg) \olg) \Y t 




3/14 



(24) 



Here Y t = h 2 is the square of the top-quark Yukawa coupling, Yf (~ 1.1) is 
the weak scale value of Y t corresponding to 'infinite' value at the GUT scale, 
and Y T is the weak scale value of the r Yukawa coupling-squared. 
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The predicted low energy values for \V C b\ are plotted in Figures 1 and 
2 (corresponding to rjp = ±1) as a function of the mass ratio \m s /md\- 
We have used the three-loop QCD and 1 loop QED beta and gamma func- 
tions to evolve the light quark masses from low energies to the top-quark 
threshold (which is also assumed to be the SUSY threshold). These run- 
ning factors, corresponding to a s (M z ) = 0.12 are (rj u , -q d , s , r) c , r) b , r) e>tl , r) T ) = 
(0.401,0.404,0.460,0.646,0.982,0.984), where the low energy scale is taken 
to be 1 GeV for (u, d, s), 1.27 GeV for c, and 4.25 GeV for b. In the square 
root in Eq. (21), only the relative negative sign (corresponding to r 2 3 > 0) 
leads to an acceptable value of \V c b\. This relative minus sign has been used 
in Figures 1 and 2. Clearly, both signs of rjp result in a consistent predic- 
tion for \V c b\. A 'central' value of V c b = 0.04 corresponds to the prediction 
m s /m d = (18, 19.5, 24) for r] P = +1 and m s /m d = (22, 21, 17.5) for r} P = -1, 
with mf hys = (165, 174, 190) GeV respectively. 

Finally, the remaining two mixing angle predictions of the model are 



\V°\ ~ 



IK° 6 I = 




1/2 \ 












e ia \ 


- f] P e 


V^o 






6 - ri2 


1/2 


6 - 


\z + r 23 


1/2 


6 + r 12 




6- 


2 Z 


- r 2 3 





(25) 



In the expression for \V® b \, the second term is numerically smaller than the 
first term, although not negligible. Taking the phase a ~ tt/2, we find that 
\V u b\ — (0.002 — 0.003). This and the \V US \ prediction resemble those from 
the Fritzsch ansatz, and are quite consistent with the present data. The CP 
parameter turns out to be J (2 — 3) x 10~ 5 . 



4. Conclusions 

We have presented a scheme for fermion masses inspired by supersymmetric 
5*0(10) in which the gauge hierarchy is implemented without fine tuning. 
This typically calls for a non-minimal Higgs sector which we exploit in de- 
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riving expression (9) for the mass matrices. Assuming spontaneous CP vio- 
lation, we are led to 6 predictions in the flavor sector which work very well, 
especially with a heavy top quark. 

There exist several ansatzes for the fermion mass matrices in the literature 
[16-21]. The Fritzsch ansatz is one of the simplest and therefore attractive. 
However, it would appear [13] that this ansatz is excluded by the recent 
Fermilab data on the top quark mass. The Georgi-Jarlskog mass matrices 
have generated renewed attention [16-17]. They lead to 6 predictions in 
the flavor sector. Owing to the up-down asymmetry in this scheme, the 
derivation of these mass matrices from an underlying theory is somewhat 
nontrivial [17]. In comparison, the scheme presented here can be obtained 
from SUSY GUTs without too much effort. There are some approaches with 
more than six flavor predictions [19], and it would be interesting to see the 
realization of these matrices from an underlying GUT or related symmetries. 

As for the neutral sector, small neutrino masses can be easily accom- 
modated in our scheme. The Dirac neutrino mass matrix is completely de- 
termined in the model (see Eq. (6)). As for the right handed (v R ) Ma- 
jorana mass matrix, we find that there are three choices which lead to a 
predictive neutrino spectrum [4,21]. They correspond to the non-zero en- 
tries in the Majorana matrix being {(11), (23), (32)}, or {(22), (13), (31)} or 
{(33), (12), (21)}, all of which result in a non-singular v R matrix. In each 
case there is a one-parameter family of solutions for the neutrino mass ra- 
tios and the mixing angles. We plan to discuss the detailed phenomenology 
of such spectra and their implications for neutrino oscillations in a separate 
paper. 
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Figure Captions 

Fig. 1. A plot of \V c b\ versus \m s /md\ corresponding to rjp = +1. a s (Mz) = 
0.12 has been used along with m c (m c ) = 1.27 GeV. The three curves 
correspond to three different values of m]? hys . 

Fig. 2. Same plot as in Fig. 1, but for rjp — — 1. 
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